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Abstract 

We study, by means of a variational method, the stability of a condensate in 
a magnetically trapped atomic Bose gas with a negative scattering length and 
find that the condensate is unstable in general. However, for temperatures 
sufficiently close to the critical temperature the condensate turns out to be 
metastable. For that case we determine in the usual WKB approximation the 
decay rate of the condensate due to macroscopic quantum fluctuations. When 
appropriate, we also calculate the decay rate due to thermal fluctuations. An 
important feature of our approach is that (nonsingular) phase fluctuations of 
the condensate are taken into account exactly. 

PACS numbers: 67.40.-w, 32.80.Pj, 42.50.Vk 
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I. INTRODUCTION 



The observation of Bose-Einstein condensation in dilute atomic 87 Rb ffl, 7 Li [[J, and 
23 Na [[J vapors last year, has created a great deal of excitement in the atomic physics 
community. Although it was speculated upon already for some time, the actual observation 
of the condensation phenomenon in these alkali gases nevertheless came as somewhat of a 
surprise, because in the spring of 1995 experiments with a trapped atomic hydrogen gas still 
held the record in achieving the necessary conditions for Bose-Einstein condensation 
The main reasons for this unexpected turn of events appear to be that alkali atoms are much 
more easy to detect and that experiments with alkali vapors can be performed at room 
temperature in contrast with the cryogenic environment that is required for experiments 
with an atomic hydrogen gas. 

The report of Bose-Einstein condensation in atomic 7 Li was also surprising for another, 
more fundamental, reason. It had namely been established experimentally that the 
effective interaction between two 7 Li atoms is attractive, or more precisely that the s-wave 
scattering length a is negative. For a homogeneous gas this implies that Bose-Einstein 
condensation cannot take place in the mechanically (meta) stable region of the phase diagram 
and is preempted by a first-order gas-liquid or gas- solid transition [[7|. Although first Hulet 
H and subsequently also Ruprecht et al. |Sj had suggested that this conclusion holds in 
an inhomogeneous situation only if the number of atoms is sufficiently large, the number 
of atoms used in the experiment indeed appeared to be too large by about two orders of 
magnitude. 

At present, this unsatisfactory state of affairs still exists, i.e. experiment claims to observe 
a condensate in a gas with a negative scattering length whereas theory seems to predict that 
this should not be possible. In an attempt to bridge at least part of the gap between theory 
and experiment, we study in Sec. |T| for zero and subsequently also for nonzero temperatures, 
the stability of a condensate of 7 Li atoms in an isotropic harmonic oscillator potential. In 
the process of this analysis, we show that there exists an interesting analogy between the 
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quantum dynamics of the condensate and the quantum mechanics of a particle in an unstable 
potential. Therefore, to make the paper more selfcontained and to bring out this analogy 
most clearly, we present first in Sec. |I| a brief summary of how the WKB approximation to 
the tunneling rate of a particle is derived by path-integral methods. We end in Sec. [TV] with 
some conclusions. 



II. TUNNELING OF A PARTICLE 

The tunneling rate r of a particle with mass m* out of a metastable minimum of a 
potential V(q) can be calculated by means of the relation 

2 2k T 

T = --Im(E ) = lim -£-Im(lnZ) , (1) 

where lm(Eo) is the imaginary part of the (analytically continued) groundstate energy in 
the metastable minimum of the potential, and Z = Tr[e~^ H \ is the partition function with 
(3 = 1/ksT and H(p, q) = p 2 /2m* + V(q) the usual Hamilton operator for the particle. This 
is a convenient starting point for our discussion, because the partition function can in the 
usual way be represented as a functional integral over the functions q(r) and p{r), i.e. 

Z = J 4q \Jm exp{-I/^ 2 *(- ip | + //( P , 5 ))} (2) 

with the periodic boundary condition q(—h/3/2) = q(hf3/2) on the coordinate but no restric- 
tions on the momentum. For the hamiltonian of interest the integral over the momentum 
p(t) is just a gaussian that can be easily carried out. As a result the partition function is in 
this case also equal to a path integral / d[q] exp {— S[q]/h} over all periodic paths q{r) and 
with a (Euclidian) action given by 

S[q]= r /2 dr(lm*(^) 2 + V(q)\ . (3) 



-h/3/2 \2 \dr 

Our task is therefore to evaluate this path integral. In general this cannot be done ex- 
actly, but one can obtain a good (semiclassical) approximation by noting that the dominant 
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contributions to the path integral are from paths that minimize the action S[q}. Such paths 
are solutions to the Euler-Lagrange equation 

m dr 2 ~ dq ' [V 

which has the same form as the classical equation of motion for a particle with mass m* in a 
potential —V(q). A periodic solution to this equation is therefore q{r) = q , where q is the 
position of the metastable minimum and obeys dV(qo)/dq = 0. Writing q{r) = q + q'{r) 
and expanding the action up to quadratic order in the fluctuations q'(r), we obtain first of 
all that 

*^/«-{4C*(H£) ,+ ^)}- (5 > 

Introducing m*ujQ = d 2 V(qo)/dq 2 , we notice that the path integral in the right-hand side is 
just equal to the partition function for a harmonic oscillator and therefore that 



[1 Ml* , ,{ m* d 2 m*ul\ , 



Afldet 



m* ( d 2 r} 
1 



-1/2 

= exp 



-^-ln(l-e-^)l j (6) 



2nh \ dr 2 

denoting the usual normalization factor due to the measure in the path integral by M . Taking 
now the limit T | (or (3 — ► oo) we thus find in this approximation for the groundstate 
energy E ~ V(qo) + Uujq/2 and for the tunneling rate T = 0. This is clearly a reasonable 
first-order result that is due to the fact that we can in first instance always approximate the 
full potential V(q) by the harmonic oscillator potential V(qo) + m*ujQ{q — qo) 2 /2. 

To obtain a nonzero value for the tunneling rate we must realize that there is another 
periodic solution to Eq. (f|) that has an action which (in the zero temperature limit) is 
only slightly different from S[qo] and therefore also gives an important contribution to the 
partition function. Using the classical analogy, this so-called 'bounce' solution ?&(t) has 
the property that the particle spends a very long time around go but in a relatively short 
time oscillates once in the potential minimum of —V(q), i.e. it bounces from q Q to q\ and 
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back to go, where q± obeys V(qi) = V(qo). In particular, for r — > ±00 it behaves as 
Qb( T ) ~ % T (^o/ a; o) e_£Jo '' r ' with t> determined by the details of the potential. 

Proceeding as before by writing q(r) = <7&(t) + q'{r) and expanding the action up to 
quadratic order in q'(r), we now find that near zero temperature the partition function 
equals 

det(-d 2 /dr 2 + lu 2 ) l 1/2 



,f3(V(q )+hu /2) 1 _|_ 



e -(S[q b ]-S[q ])/h I ? ^ 



_det(-d 2 /dr 2 + io b (r) 2 ) 

introducing the quantity u;&(t) by means of 7??*u;&(t) 2 = d 2 V(qb(r))/dq 2 . Adding in a similar 
manner also the contributions from paths with an arbitrary number of bounces, we find that 
the series exponentiates and hence that 



E ~ V(q ) + ^ - hm k B T 



det(-d 2 /dr 2 +wg 



11/2 



e -(Sfe]-S[9o])/» 



_det(-rf 2 /rfr 2 + ^(r) 2 ) 

We expect the third term in the right-hand side to represent the tunneling rate out of 
the metastable minimum and therefore to be purely imaginary. This expectation is indeed 
correct, because the operator —d 2 /dT 2 + uJf J (T) 2 turns out to have a negative eigenvalue. This 
is most easily understood from the fact that Eq. (f|) shows that dqb(r)/dT is an eigenfunction 
of this operator with an eigenvalue equal to zero. Since this eigenfunction has one node, we 
know from our experience with the Schrodinger equation that there must be an eigenfunction 
without nodes that has a lower, and therefore, negative eigenvalue. 

However, the presence of an eigenvalue equal to zero appears to give an infinite result 
for the ratio of determinants in Eq. (|8]). Fortunately, this is due to an improper treatment 
of the zero mode in the calculation of the path integral over the fluctuations q'{r). Since 
qb(r — t ) = qbij) — Todqb(r)/dT + 0(tq ), we note that this zero mode is just associated with 
a translation of the 'bounce' solution (#,(t) and hence that the square root of the ratio of 
determinants must be proportional to the total time interval Hf3 = h/hsT. A detailed and 
beautiful analysis by Duru et al. actually shows that 



Combining this with Eq. @ we finally arrive at 

r = -hm(Eo) = JH^L e -(SM-sM)/n (1Q) 

where 

is recognized as the usual WKB expression for the exponent of the tunneling rate. 

III. TUNNELING OF A CONDENSATE 

We now turn to the problem of the stability of the condensate in a magnetically trapped 
gas of 7 Li atoms. We consider here only the experimentally relevant case of a large number 
of particles iV ^> 1 in a large trap. Quantitatively, the latter means that ry/£ <C 1, where 
ry is the range of the interatomic interaction, t = ^nh/mu! is the spatial extent of the 
one-particle ground state in a harmonic oscillator potential with level spacing hu, and m 
is the mass of the 7 Li atoms. (In the case of 7 Li the above condition also implies that 
\a\/£ <C 1.) Moreover, we will always require that the density n in the center in the trap is 
such that the gas parameter nry <C 1. This basically leads to an upper bound on the total 
number of particles that is always satisfied by the experiment of interest ||. 

A. The case T = 

At zero temperature and in an external trapping potential V ext (x) = mu) 2 x 2 /2, the 
effective hamiltonian for the condensate wavefunction ip(x,t) and its canonical momentum 
7r(x,t) = ihtp*(x,t) is given by 

F[tt,V]= / dx r (-—V 2 + V ext + T2B ^ ° ; 0) M 2 ^) i) (12) 
J \ 2m 2 J 

since then the Hamilton equations exactly reproduce the nonlinear Schrodinger equations 
for ip(x,t) and i/j*(x,t), respectively. For example, we correctly find that 
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^(-^ 2 + V^\x)+T 2B (6,6;0M(x,t)A^(x,t) , (13) 



with T 2B (0, 0; 0) = Aitafi 2 /m the effective interaction between the atoms. It is important to 



mention here that, due to infrared divergences in the theory of the dilute Bose gas [11], this 
nonlinear Schrodinger equation is only valid if the energy cut-off provided by the trapping 
potential is sufficiently large. More precisely we must require that hcu > Ann\a\h 2 jm. Fortu- 
nately, this correponds precisely to the conditions under which the condensate is metastable 
as we will find out shortly. 

After this brief discussion of the condensate wavefunction as a classical field, we can now 
turn to the quantum fluctuations of the condensate. As in Sec. |TT], we must then consider 
the partition function Z = Trie -1311 ] at zero temperature. This function can be written as 
the functional integral 

Z = j '(#*]<#] expj-iyV (j dx i)*h^ + H[ih^\i)]j^ (14) 

over the periodic fields ip(x, r) and ip*(x, t), which is the direct analog of Eq. (0). To proceed 
as in the case of the quantum mechanics of a particle, we should now integrate out the 
momentum field. This is, however, not helpful for our purposes because the hamiltonian does 
not have a momentum independent part. Consequently, it seems that a stability analysis of 
the condensate cannot be performed in the same way as in Sec. ||. 

The way out of this dilemma is found by noting that the instability of the condensate 
is, just as in the homogeneous case, a result of density fluctuations that lead to a lower 
energy because of the effectively attractive interaction between the 7 Li atoms. It is therefore 
advantageous to perform a canonical transformation by means of the relation ip = J~pe %x 
and to use the density field p(x, r) and the phase field x(af, r) to calculate the partition 
function. Indeed, a simple substitution together with the periodicity of p(x, r) shows that 
the partition function is equal to the functional integral J d[p] J d[x] exp{— S[p, x\/Ti} with 
an action 



that is quadratic in the phase field r). Hence, we can now immediately integrate over 
this field. 

There is an important point to be made about this integration, which reflects the fact 
that if the original fields ip(x,r) and ip*(x,r) are periodic, the phase field r) is only 
periodic up to a multiple of 2n. To calculate the partition function correctly, we must 
therefore first integrate over all the fields i~) obeying the boundary condition oo) = 
x(x, — oo) + 2nj and subsequently sum over all possible integers j. Clearly, this change in 
boundary conditions affects only the zero- momentum part of r). As a result we have 
to consider the sum 



dr 
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first, where we made use of the fact that the total number of particles in the condensate 
No(t) = J dxp(x, r) and the boundary condition Xo(°°) = Xo(~ 00 ) + ^3 011 t ne integration 
is implicitly assumed for each term in the sum. Performing a partial integration on the 
integral in the exponent, this sum becomes equal to 

r v dT - 

because the integration over Xo( r ) then simply leads to the constraint of a constant number 
of particles in the condensate, i.e. N (t) = N . At zero temperature and under the conditions 
that the nonlinear Schrodinger equation is valid, this actually implies that N (t) is equal to 
the total number of particles N. In addition, the sum J2j e 2mNoj equals J2j — j) and 
thus requires that N is an integer. In this manner we see explicitly that the integration 
over the zero-momentum part of p(x, r) is effectively only a sum over the total number of 
particles and, most important for our purposes, that the partition function at a constant 
number of particles is given by the functional integral / d[p] j d[x] exp{— S[p, x]/h} over a U 
nonzero momentum components of the density and phase fields. 



The integration over the nonzero momentum part of x(x, r) is readily accomplished by 
introducing the Green's function for the phase fluctuations G(x, x 1 ; p), that obeys 

n 



m 



((Vp) ■ V + pV 2 ) G(x, x- p) = 5{x - x') . (16) 



In terms of this Green's function the formal solution to the Euler-Lagrange equation 
S3\p,x]/Sx = reads 

X{x,r) = -i J dx 1 G(x,x l ;p) 9p ^ , (17) 



which after a substitution in Eq. (|T5|) gives the following effective action for the density field 

h dp(x,r) dp(x' ,r)\ 



S[p] = JdrJdxJ dx' f-| P[ * ,T) G(x,x';p)- 



dr 



+ f iT f iS (JL iVpr + v ^ p+ ^M>A. (18) 



imp 2 

This is the desired analog of Eq. (|3]) that we will now use to determine the stability of the 
condensate at zero temperature. 

Our first task is to see whether a metastable condensate is possible at all in a trap. This 
question can be answered by considering only the 'potential energy' 

vw = /«f^ + n + ^) d9) 



imp 2 

for time-independent density configurations. Because we are unable to analytically consider 
all configurations, we proceed in a variational way [12|] and investigate here only the gaussian 
profiles 

p(f; g) = AT 1 exp(-|J) . (20) 

The physical reason behind this choice is that we expect the shape of the condensate, if it 
is metastable, to be close to the square of the one-particle groundstate wavefunction of the 
trapping potential and, therefore, to be reasonably accurately described by a gaussian. We 
will see shortly that this expectation is indeed correct. 
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Substituting the above gaussian into the expression for the 'potential energy', we find 
that 



Hence, for \a\ = the potential V(q) has an absolute minimum at q — = £/^/2n 

and the 'potential energy' at this minimum is 3Nhu/2. Clearly, this corresponds exactly 
to an ideal Bose condensate. For \a\ 7^ the 'potential energy' is always unbounded from 
below because V(q) — > —00 if q J. 0. This implies that the condensate is in general unstable 
and tends to collapse to the density profile lim ? | p(x; q) = NS(x). Of course, in a realistic 
system the collapse to this density profile is ultimately prevented by the hard core of the 
interatomic interaction as the density reaches a value of the order of 0(1 fry). Although 
our theory certainly breaks down at these densities, the approximation NS(x) for the final 
density profile is nevertheless rather accurate for the lengthscales of interest, due to the 
condition ry/£<^ 1. 

The most important feature of Eq. fl2"ip is, however, that if the condition 



is fulfilled, the potential V(q) has a metastable minimum. This result can directly be com- 
pared with the work of Ruprecht et ai, who find by a numerical integration of the nonlinear 
Schrodinger equation Eq. ( |T3"D that a stable solution can only be obtained if N\a\/£ < 0.23 
M. Our variational calculation, therefore, gives an upper bound that is only 16% too high. 
Apparently this is the amount of accuracy that one can obtain for the condensate energy by 
considering only the density profiles p(x;q). 

Having arrived at the conclusion that a metastable condensate is possible at zero tem- 
perature if the number of particles in the gas is sufficiently small, we now want to calculate 
the lifetime of the condensate due to macroscopic quantum fluctuations. This involves also 
the evaluation of the 'kinetic energy' 




(21) 





(23) 
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Again an exact treatment of this part of the effective action is unfeasible. However, our 
previous discussion of the 'potential energy' suggests a reasonably accurate approximation 
that amounts to the promotion of the variational parameter q to a real dynamical variable 
q(r). Hence, we assume that for a description of the tunneling process the most important 
configurations of the condensate are given by p{x, r) = p(x;g(r)). The reason for this 
assumption is threefold. First, we have already seen that the profile p(x; q) gives an accurate 
description of the metastable minimum of the potential V[p\. Second, one expects that the 
minimum energy barrier for the tunneling process is associated with a rotationally symmetric 



configuration [T3| that, because of the diluteness of the gas, can be represented by a gaussian 
with the same amount of accuracy as the metastable minimum. Third, the profile p(x; q) 
describes in the limit q { also the completely collapsed state of the system. In combination, 
the configurations p{x\q(r)) thus seem to provide a reasonable interpolation between the 
initial and final density profiles of the tunneling process. 

Evaluating the 'kinetic energy' for the variational ansatz p(x, r) = p(x; q(r)) and adding 
Eq. ( |2"TD we easily find that the quantum dynamics of the condensate is determined by the 
action 

... , ) 2 

which differs only by an overall factor iV from the zero temperature limit of Eq. (JJ). The 
effective mass m* is difficult to calculation exactly due to the inhomogeneity of the system, 
but can be estimated by noting that the dominant contributions to the 'kinetic energy' come 
from the region near the center of the trap. In this region the Green's function for the phase 
fluctuations is well approximated by 

Tfi I 

G(x, x 1 ; p) ~ -V- - - (25) 

Mp(0;q)\x-x'\ v ; 



S[q]=N I dr\ l -m*(^\ +V(q)\ , (24) 



and we obtain 



* m 
m ~ — 



^5/2 Jdyjdy' (j - j?) exp {-?)-!— (j - y>^ exp (-^ 2 ) ~ 0.27m . (26) 
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Although the exact value of the effective mass is presumably somewhat larger than this 
result, it is in any case of the order of the atomic mass m. 

Applying now the final result of Sec. |I| (replacing only H by fr/N) we immediately find 
that the decay rate of the condensate is equal to 

r = exp {_™ £ dq j 2m , mq) m 

and typically very small for a number of condensate particles N ^> 1 but still sufficiently 
small to fulfill the condition of metastability. Note that both requirements do not exclude 
each other, because we are only considering traps for which \a\/£ <C 1. 



B. The case T^O 

We have seen that at zero temperature a metastable condensate is possible if the number 
of 7 Li atoms is less than N m = 0.23£/\a\ ~ 1440. (The scattering length is — 27ao and 
the effective isotropic level splitting ^o^o^u^) 1 / 3 for the trap is 7nK 0.) In the Rice 
experiment, however, one reports the observation of a condensate for temperatures as high 
as AOOnK and total numbers of atoms as great as 2 • 10 5 . Clearly, for these temperatures 
the previous discussion is inadequate and needs to be modified before we can draw any 
conclusions about a possible discrepancy between theory and experiment. 

The generalization to nonzero temperatures is essentially only useful in the case that the 
number of particles iV ^> 100, since then the critical temperature T c is almost equal to the 
critical temperature 

of the ideal Bose gas and obeys ksT c ^> tkv so that we are close to the thermodynamic limit. 
As a result also the thermal wavelength A c = \j2nTi 2 /mk,BT c <C t. To understand why 
this is important, we must realize that at nonzero temperatures the nonlinear Schrodinger 
equation for the condensate wavefunction is 
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ih^(x,t)= (-|V + ^(f)+T 2B (0,^ , (29) 

with n'(x,t) the local density of atoms which are not in the condensate ||14]| , and that the 
total number of noncondensed particles is of 0((/csT '/huj) 3 ) because n'(0, t) is of 0(1/A 3 ) 
and the size of the noncondensed cloud is of 0(£ 2 /A). Therefore, for temperatures that 
do not obey ksT ^> hu>, the term 2n'(x,t) in Eq. (p9|) is negligible compared to \ip(x, t)\ 2 . 
For temperatures T ^> huj/ks this does not have to be the case. However, under these 
conditions we can replace 2n'(x, t) by a constant since the size of the noncondensed cloud is 
now much larger than the size of the condensate. Performing then the same analysis as in 
Sec. Ill A , we find that the potential V(q) is just shifted by a constant. 

We thus conclude that for A" 3> 100 a metastable condensate can exists as long as 
N < N m , or sufficiently close to the critical temperature. Moreover, the decay rate of the 
condensate can at all temperatures be calculated by means of the effective action 



S M= N °lZ dT \2 m '(£) 2 + V{q) ) ' <30) 



Hence, for T <C Tooo/hs quantum fluctuations dominate and the decay rate is given by 
Eq. (p7D with A" replaced by Ao- For T ^> huo/kB, which is the relevant temperature interval 
if N ^> N m , the decay of the condensate is dominated by thermal (or classical) fluctuations 
and the decay rate is proportional to the Boltzmann factor e" N °^ v ^ qm ^~ v ^'' kBT , where q m 
denotes the position at which V(q) has a maximum. The prefactor of the exponent is 
always a difficult point, but can be estimated by noting that in the classical limit the action 
of Eq. (|3TiD implies a probability distribution 

P(p, q) ^ exp | -p + N (V(q) - Vfa,))) } (31) 

for the momentum p and the coordinate q of a fictitious classical particle. Using this prob- 
ability distribution we then easily find from the average flux over the energy barrier that 

ro = ^exp{-^(K(g m )-^o))} • (32) 
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Notice that this decay rate is only small if N (V(q m ) — V(qo)) <C ksT. Therefore, a long- 
lived metastable condensate is only possible sufficiently close to the critical temperature, if 
the total number of atoms obeys N <C (N m ) 3 . 



IV. CONCLUSIONS 

In view of the results obtained above it seems that we have arrived at an explanation of 
the experiment with atomic 7 Li, since we have shown that a long-lived metastable condensate 
is possible if 1 <C iV < N m . Unfortunately, this is only partly true because we have not yet 
considered the mechanical stability of the noncondensed cloud. In particular, we have not 
shown that the gas is stable (or metastable) at the critical temperature and will, in contrast 
to what occurs in the homogeneous case 0, not immediately collapse to a dense liquid or 
solid phase. In our opinion, this is at present the most important question which remains 
to be answered before we can speak of a theoretical understanding of the Rice experiment. 

Finally, we need to discuss the important observation by Dalfovo and Stringari that 
a condensate with a vortex line through the center of the trap has a value of N m that is 
significantly larger than 0.23£/|a| , due to the fact that the local density of the gas must vanish 
at the point where the phase singularity occurs |L5]. In principle, we can also consider the 



decay of a condensate with a vortex in the framework of the theory presented in Sec. [TTTj. We 
should then first put x(af, r) = Xv(%, T ) + x'i^i T )> where Xv(x, T ) is the phase configuration 
of a vortex, before we integrate over the nonsingular phase fluctuations x!(x,t). The main 
difficulty, however, will be to find an accurate solution of the coupled dynamical equations 
for the density profile and the position of the vortex. Although such a calculation would 
certainly be interesting in its own right, it clearly does not invalidate the main conclusion 
of our paper that the most important unresolved issue is the stability of the gas without a 
condensate. 



After completion of the work presented in Sec. [Ill A| of this paper, we received a preprint 



by Kagan, Shlyapnikov, and Walraven ]IB[ in which these authors also consider the stability 
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of an inhomogeneous condensate at T = 0, and arrive at qualitatively the same results. 
However, their calculation of the decay rate is completely different from ours and, in partic- 
ular, does not take phase fluctuations of the condensate into account. As a result, it appears 
not to correspond to a WKB treatment of the tunneling process. 
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